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Abstract — In this paper, by way of three examples - a fourth 
order low pass active RC filter, a rudimentary BJT amplifier, 
and an LC ladder - we show, how the algebraic capabilities of 
modern computer algebra systems can, or in the last example, 
might be brought to use in the task of designing analog circuits. 
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I. Introduction 

For those about to embark on the study of this paper, 
let us start with a warning: This article is not a report on 
ground-breaking results. Neither will we present spanking 
new circuits, ripe for patent. Rather in spirit of SM^ACD's 
explicit policy to be "a forum for the discussion of new 
ideas and methodologies", we would like to advocate the 
(re)introductior(3 of computer algebra into the design process 
of analog circuits. In order to do this, we present some 
examples, where computer algebra (together with numerics) 
was used in the design, especially the sizing of the given 
circuits. 

These examples do not represent a systematic approach, but 
they should be considered as tentative steps towards getting 
to know the capabilities (and some of the shortcomings) of 
computer algebra systems (CAS), and making these tools bring 
forth fruit in the daily work of electrical engineers. 

II. The general ansatz 

One of the principal problems in analog circuit desigrH is 
to determine element values in a given circuit in such a 
way that the resulting network function is in accordance 
with some previously given performance specifications. In 

This aiticle first appeared in: SM^ACD'08 Proceedings of the Xth In- 
ternational Workshop on Symholic and Numerical Methods, Modeling and 
Applications to Circuit Design, Erfurt, Germany, October 07 to 08, 2008; 
pp. 127-134. Due to the low distiibution of these proceedings, the author has 
decided to make the article available to a larger audience through the arXiv. 

The author's current (as of December 31st, 2010) address is BismarckstraBe 
20, D-38518 Gifliorn, Germany. Cun'ent e-mail: e.gerbracht@web.de 

' Reintroduction, because some symbolic analysis tools were and still are 
closely associated to existing computer algebra systems, although nowadays 
they do not make much use of the advanced algebraic capabilities (it is 
debatable if they ever did), but rely more on their numeric proficiencies. 

2 See e.g. (T]. 



more mathematical terms: suppose the network function is 
demanded to be a rational function of form 

ao + flis + • • • + a^s™ 



fdes{s) 



(1) 



K + his^ h6„s" 

with the numbers m, ai, . . . , Om, bi, . . . ,bn determined by 
the initial specifications. 

Then we need to find an analog circuit F built from 
elements, which themselves are described by parameters 
Pi, . ■ . ,Pk, such that its network function 



Hr{s) 



(2) 



(4) 



E"=o B^{pi,...,pk)s'' 

where the terms and are polynomials in pi, . . . ,pk, 
satisfies 

Hr{s) = fdesis). (3) 

In other words: at the end of the day, many design problems 
result in a system of nonlinear equations 

Af,{pi,...,Pk)-af,= 0, 
B^{pi,...,Pk) -b^ ^ 0, 

1 < < m, and 1 < ly < n for variables pj, . . . ,pk- 
If indeed the terms A^ , Bi, are polynomials^ with coefficients 
in some number field, the set of solutions, otherwise also 
known as the zero set of (|4]l, forms what is called an affine 
algebraic variety, an object which belongs to the mathematical 
subdiscipline of Algebraic Geometry. Thus when trying to 
solve these equations, we could profit both from the centuries 
of experience, mathematicians have gained in this particular 
field, and from the more recent results and algorithmic instru- 
ments developed - the concept of Grobner bases, as well as the 
Buchberger algorithm and its "relatives" which are realized in 
nearly all of the current computer algebra systems. The result 
would be exact ("symbolic") descriptions of the sought after 
solution sets, from which we could deduce numeric results 
afterwards. 

Now, when we take a look at the advances in symbolic 
methods applied to circuit design, we can make a peculiar ob- 
servation: nowadays it has become nearly routine to calculate 
any given network function as a multivariable rational function 
of the complex frequency s and the symbolic circuit elements 
(even though there is still the problem of the combinatorial 

^Thus for example we do not talk about the problem of biasing a transistor, 
since then we would have to solve equations containing the exponential 
function (as in the Ebers-Moll model of the BJT), or equations which contain 
roots, and moreover are piecewisely defined (as for FETs). 
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explosion for large circuits). However, when it comes to 
finding solutions for the nonlinear equations arising from the 
task of sizing a particular circuit, even if it is relatively small, 
researchers, as well as designers, fall back upon numerical 
means and various approximative methods with their various 
drawbacks. 

III. An Example - Sizing an active 4th order low 

PASS RC filter 
In section 6.3.3. of ^ the following filter circuit was 
presented: 



Fig. 1. Circuit topology for fourth order low pass RC filter 

The aim there was to construct a fourth order Butterworth filter 
with this topology. 

Using one of the available symbolic analysis programfl we 
see that the circuits symbolic transfer function is given by 



(5) 



with il'(s) = l + ais + a2s'^ +a3S^ + a4s'^ being a polynomial, 
and 

ai = {(l-K){Ci+C2)+C3+Ci)R4. + i{l-K)C2+C3+Ci)R3 + 
{{l-K)C2+Ci)R2+CiRu 

a2 = C2CiRiR2 + {C3+C2)CiRiR3 + {Ci+C3+C2)CiRiRi + 

C3((l-K)C2+Ct)R2R3 + ((Ci+C3)((l-K)C2+Ct)R2Ri + 
Ci{C3 + il-K)C2+Ci)R3R4, (6) 

a-S = C3C2ClRlR2(R3 + Ri) + CAC2ClRl{R2+R.3)Ri + 
C4 C3 Cl ( fll +fl2 ) ^3 fl4 + ( 1 - -ff ) C4 C3 C2 fl2 -R.3 -R4 , 

04 = CiC3C2CiRiR2R3Ri- 

In order to specify element values in such a way that H{s) be- 
comes the transfer function of a fourth order Butterworth filter, 
we can proceed from these data by matching the coefficients 
with those numerical values, which can be found in any table 
for filter design, and then try to solve the resulting equations by 
numerical means. Since we only have four equations for nine 
unknowns, there is the need to reduce the existing degrees of 
freedom by more or less arbitrarily attaching numerical values 
to some of the element values. In fact, that is what was done 
in 121, where K was set equal to 2 and all the Ci were set 
to 1 (all element values in this example are considered to be 
normalized); afterwards some variant of the Newton-Raphson 
method was used, to deduce one single solution. 
The disadvantages of such an approach are obvious: Using 
rounded values (from precalculated tables) introduces errors 
right from the start. Both, the assignment of (quite arbitrary) 

''For the analysis of this and all other circuits within this paper, we use 
SapWin 3.0, because of its affordabiLity and its easy availability through the 
world wide web |3j. 



numbers to arbitrarily chosen parameters, and the choice of a 
starting vector for the Newton-Raphson method might hinder 
from finding a solution, or even worse, might make a solution 
impossiblqj 

Furthermore, questions like "Do alternatives exist? How 
many are there?", cannot be answered either, and changing 
any of the numerical parameters forces us to redo most of the 
calculations. 

When we tried to solve this design task, first of all we 
started with a slightly different approach, using the Feldtkeller 
equation 



D{s) ■ D{-s) = 1 



(7) 



for the special case of the fourth order Butterworth filter Thus 
we got rid of the precalculated approximate values for the 
Oi, by using instead the defining algebraic equations given by 
matching the coefficients from 

D{s) ■ D{~s) ^ 1 + (2a2 - ai)s^ + {aj - 2aia3 + 204)5" 
+ (20204 - 03)5'^ + 045^ 

with 1 + s^. Consequently we enlarged the set of equations 
^ by adding 

2o2 — o^ = 0, O2 — 201O3 + 2o4 = 0, 202O4 — 03 = 0, 04 = 1 

(8) 

to it, and regarding oi, . . . , 04 as further unknowns. 

Even though in this paper we propose the usage of the 
algebraic capabilities of a CAS, at this point we used Math- 
ematica (Version 6.0.1.0) as a numerical engine - not in 
order to approximate one solution, but to find approximates 
for all solutions. Moreover as was done in we also fix 
K = 2, and Ci — 1 for 1 < i < 4. Since obviously 
the enlargement of the set of equations results in a bigger 
set of solutions, we were left with the final task of "post- 
processing" this zero set, by singling out those particular 
solutions for which oi, . . . , 04 determine a polynomial D{s), 
all of whose zeroes have negative real part (= filter stability), 
and such that all the Ri are non-negative (= realizability 
by classic elements). Indeed, we observe that from all the 
possible solutions calculated, only nine satisfy the restrictions 
concerning the o^, and again only one of these, 

K = 2, Ci = C2 = C3 = C4 = 1, 

Ri w 0.133933818297194652631087580090, 

R2 ~ 3.893036697318392402871746149006, 

R3 w 2.479192111455558403082198766696, 

Ri w 0.773590398536329977043175927841. 

results in the Ri all being positive. 

So, where does computer algebra enter this picture? 

First of all Mathematica's NSolve-procedure, which we 
used in the above task, at its core heavily depends on algebraic 
algorithms: the fact that it is able to calculate all solutions 
stems from the algebraic variety, which was defined by ^ and 
(Ull, being zero dimensional. That is, it consists of finitely many 

^E.g. this is the case, when we use the preassigned values from |2|, and try 
to find element values, such that the above circuit attains the characteristics 
of a fourth order Chebyshev filter. 
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points. Again we can be sure of this fact, only because there 
is an algorithm - the Buchberger algorithr^ - which is used 
by Mathematica as a subprocedure of NSolve and which for 
any system of algebraic equations produces particularly well- 
behaved equivalent systems of equations - so-called Grobner 
bases - which allow much more insight into the properties of 
the zero set, and may be used for a subsequent elimination of 
variable^ 

Thus e.g., when we apply Mathematica's 
GroebnerBasis-command to Q and (O, with the 
values of K,Ci, . . . , C4 fixed as above, and the variables 
ordered lexicographically by ai < ... 04 < Ri < ■ ■ ■ < R4 
(which amounts to saying that first ai is to be eliminated, 
then 02, etc.), then the resulting Grobner basis contains the 
polynomial 



P 



1- 160i?|4 
-920274816 i?f- 



13872 i?^ 



788512 + 31505120 i?f 



20065991808 iif 



-328437088768 Rf + 4000414289152 Rf 



-35535204282368 i?f 



223781766674432 i?f 



-956822102532096 i?f -I- 2535921430958080 Rf 
-3050522934050816 i?f - 2341746368053248 i?f 



-10182414501412864 



13996296348106752 



■68 



1806331484831744 i?f 
2888816545234944 i?P 



in i?4 of degree 72, which implies that there are only 
finitely many possible values for R4. Furthermore this Grobner 
basis contains polynomials pi(i?i,i?4), ^2(^2,^4), and 
P3(i?3,i?4), each of which being linear in the respective Ri, 
and otherwise consisting only of a high degree polynomial in 
R4. Since these polynomials as part of a Grobner basis have to 
be zero, too, consequently each zero of P determines precisely 
one solution vector {Ri, . . . , R4). 

Those solutions, which result in the additional specifica- 
tions for the ai to be fulfilled, originate from a factor of 
the polynomial P of degree 18, which can be calculated 
using Mathematica's Factor-procedure, when we allow the 
coefficients to be sums of integers and integer multiples of 
\/2. This factor polynomial is 



-10 + 7V2 + (-12 + 8V2)Rl 
+ (520 - 352\/2)i?^ + (-2824 
+ (-10816 + 6656%/2)i?f 
+ (-8288 + 5952V2)i?^''- 
+10368i?F. 



- (280 
1900 



196V2)i?: 



2 i?2 



(-13904 + 7368 V2)i?f 
(1280 + 6592V2)Rl 



To finish this section, let us remark that by modifying the 
Feldtkeller equation to 

D{s)-D{s)^l + T4(s)-T4{-s), 

where T4(s) denotes the fourth order Chebyshev polynomial, 
and slightly changing the initial values for the Ci, after some 
trial-and-error we were able to calculate the following set of 
parameter values, with which we can make the initially given 
circuit into a Chebyshev filter 



K 
Ri 
R2 

i?3 
i?4 



2, = 1, C2 = 2, C3 = 2, C4 = 1, 
0.263638090854794185461593787592, 
0.624164765447879000316525786179, 
2.645185226758545312744750592839, 
3.249013399873649633437777451232. 



The question remains for each filter characteristic: which 
set of prescribed values for if, Ci , . . . , C4 imply solutions for 
i?i,...,i?4, which make sense in an electrical engineering 
point of view, i.e., for which positive JsT, Ci, . . . , C4 do we 
get positive Ri, . . . , R4I 



IV. Another Example: Sizing a BJT Amplifier with 

RESPECT TO BREAK FREQUENCIES 

From 16j - a students' textbook in analog electronics - we 
take our next example: a rudimentary BJT amplifier circuit 
with only a few additional classic elements for biasing and 
sizing purposes (see figure (llVb ). The biasing resistances are 



Fig. 2. Rudimentary topology for a BJT amplifier 

given by Rb = 310kf7, and Rc — 2kil. The biased transistor 
is modelled by concrete h-parameters 



Substituting T = R^, we see that its structure is effectively 
determined by a polynomial of degree 9, which results in the 
number of admissable solutions being nine. 



H 



All 


hi2\ 




U21 


/122/ 







35 • 10/zS 



(9) 



^The Buchberger algorithm can be seen as a generalization of both the 
Euclidean algorithm (to the multidimensional case) and the GauB algorithm 
(to polynomial equations of degree > 1). 

'in most cases, especially if the zero set is zero dimensional, there exist 
Grobner bases which are in triangular form, a property which generalizes 
the analogous property of sets of linear equations, which can be brought to 
triangular form using the GauS algorithm. For more details see (4), (5). 



and its capacitances Cbc = HpF, and Cbe = 320pF. 

With denoting its transfer function by H{s) ~ ^(fj; 
symbolic analysis with Sap Win 3.0. gives 

N{s) = Ceh2lRBRc ■ S - CBC'CehnRBRc ■ 



a full 
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and D{s) — uq + ais + a2S^ + 035^, with 

do = -/lll--RB+hi2/l21-Rc-''ll'l22-Rc— '122-Rb-Rc, 
0,1 — —CBchllRB-CBEhllRB-CJlllRB-CahllRc 

— CBcf^llRc ~CaRBRc^CBcRBRc+CBC^12RBRc 

— CBch2lRBRc+CBchl2h2lRBRc+CBEhl2h2lRBRc 
+ dhi2h2lRB Rc - C Bchllh22RBRc —C BEhllh22RB Rc 

-C^hiih22RBRc-C^hiiB.,-C,RBRs+CJii2h2iRcRs 

-C^hiih22RcR^-C^h22RBRcRs, 

02 = -CaCBchllRBRc-CaCBEhllRBRc-CBcCBEhllRBRc 

-CaC^hiiRBRc-CBcC^hiiRBRc-CBcC^hiiRBRs 
-C'BEC^hiiRBRs-CaC^hiiRcRs-CBcC^hiiRcRs 

-CaC^RBRcRs-CBcC^RBRcRs+CBcC^hi2RBRcRs 

-CBcC^h2lRBRcRs+CBcCehi2h2lRBRcRs 

+ CBEC^hi2h2lRBRcRs-CBcCehiih22RBRcRs 

-CBEC^hiih22RBRcRs, 

03 = -CaCBcC^hiiRBRcRs-CaCBEC^^hiiRBRcRB 

— CbcCbeC^HwRbRcRs- 

With regard to the concretely given transistor parameters at the 
bias point, the nominator and denominator of H{s) become 

N'{s) = (290625000000Ces(-1000000000000+ 77s)), 

and 

D\s) = 1623139843750000000 
+3590505265625 s 
+3033906250000000000000 CaS 
+1088390625000000000000 CgS 
+1623139843750000000 CeRsS 
+7161 

+673378125000000 CaS^ + 22378125000000 CeS 
+2034375000000000000000000CaCeS^ 
+3590505265625Cei?sS^ 
+3033906250000000000000CaCei?^s^ 
+7161Cei?sS^ + 673378125000000CaCei^sS^ 

respectively. Thus, e.g. setting Rs = Oil, Ca = 300pF, and 
Cfc = 0.1/iF, we get 

, 465000s(-1000000000000+77s) 

His) — , 

25970237500000000000+1813435834250000s+1015651791s^ 

(10) 

the poles being (approximateljo 

pi « -14437.758, and p2 ~ -1, 771051.964. 

This results in the Bode plot shown in figure [3] 

We will now use this example as a problem in circuit design. 
Our task shall be: with the bias point remaining unchanged, 
find parameters Rs,Ca, and Ce in such a way that a new 
transfer function H{s) is achieved, which has the same zeroes 

^Clearly, due to the nominator being a quadratic function, we would also 
be able to calculate exact values. 
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Fig. 3. Bode Plot of initially given BJT amplifier, according to (TO) 



as H{s), and possesses two (new) poles pi andp2- Introducing 
a new parameter k (since due to the simplicity of our design 
topology, we cannot specify the maximal attenuation of H, 
as we will see below), we thus are asked to determine the 
parameters k,Rs,Ca, and Ce in such a way that 

N'is) fc • s • (77s - 1000000000000) 

„ , , = His) = ; ^ -, ^ . (11) 

This is equivalent to 

iV'(s)-(s-pi)-(s-p2)-£''(s)-fc-s-(77s-1000000000000) = 0. 

(12) 

Factoring out common terms on both sides of the — 
sign, we are left with a third order polynomial in s 
with coefficients being integer polynomials in the variables 
{k, Rs, Ca, Ce,Pi,P2] ■ To solve (fTZt we need to set these co- 
efficients to zero, which consequently results in the following 
system of four algebraic equation^ (the z-th equation resulting 
from the coefficient of s'): 

= 1623139843750000000 k - 290625000000 CePiP2, 
= 3590505265625 k + 3033906250000000000000 C^fc 

+ 1088390625000000000000 C^k + 290625000000 C^Pi 

+ 290625000000 CeP2 + 1623139843750000000 Cefcii^, 
= -290625000000 Ce + 7161 /c (13) 

+ 673378125000000 Ca/c + 22378125000000 Cefc 

+ 2034375000000000000000000 CaC^k 

+ 3590505265625 CekRs 

+ 3033906250000000000000 CaC^kR^, 
= 7161 C'ekRs + 673378125000000 CaCefci?s- 

Alhough this system is nonlinear, using its procedure Solve, 
Mathematica is able to calculate four symbolic solutions for 
k,Rs,Ca, and Ce as functions of pi and p2, voicing the ex- 
plicit caveat: "Equations may not give solutions 
for all "solve" variables .". Nevertheless, indeed 
there are only four different solutions. 

'This is why we had to introduce the additional parameter k, since we 
consider pi and p2 to be fixed. 
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The first is irrelevant for the concrete problem at hand from 
an electrical engineering point of view, since it is given by 



Ce = and 0, 



(14) 



and thus results in a trivial transfer function. 

The second one, though it leads to nontrivial terms for Ce, k 
and Rs, furthermore implies 

11 



Ca 



(15) 



1034375000000' 
which again has to be considered absurd from an engineering 
point of vie w0 

Finally, there are two solutions left which lead to two 
different valid sets of design parameters, both of them having 
in common that Rs = 0. 

Here we will not present the explicit formulas for Ca,Ce,k, 
but we will give a numerical example, and we will further 
show, how the formulas for Ca (and beyond this our claim, 
that there are only four solutions) can be derived with the help 
of Grobner bases. 

First, for the numerical example let pi — —10, and p2 — 
— 1000. Then the two remaining approximate solutions of the 
above set of equations are 

(Ca, Ce, k) w (53.4989AiF, 1.49102/^F, 2.66969 • 10"^) (16) 

and 

(Ca, Ce, k) « (53.5000AiF, 149.12881mF, 2.67017 • 10"^). 

(17) 

Indeed, the corresponding Bode plots for the respective 
circuits can be seen in figure 4. 




Fig. 4. Bode plot of BJT amplifier circuit with original parameters (dotted), 
parameters 116) (dashed), and (TT) (thick) 



Finally, applying Mathematica's GroebnerBasis-comm- 
and to the right-hand sides of equations (fT3] l. with the order 
of variables chosen to be fc < i?s < Ca < Ce, results in a 
Grobner basis consisting of 11 polynomials with exactly one of 
them containing only the variables Ce,Pi,P2- This particular 
polynomial is of degree four in the variable Ce- Keeping in 
mind that the zero set of the original polynomials constituting 

'"Unless we allow for NICs (negative impedance converters) to be added 
to the initial circuit topology. 



(fT3] l is equal to the zero set of its Grobner basis, we thus see 
that there truly are only four possibilities for the values of Ce- 
Factoring this particular polynomial (again with the help of 
Mathematica) we get the following factors: 



Ce 



and 



(18) 



1182300335490970802500000000000000 
~ 2591834572818828634150000000pi 
- 2591834572818828634150000000 p2 
f 5681810493667293811609piP2 
f 1737403713997011459000000000000 pip2 Ce, 



and finally 

100852820307500000000000000 

- 22384370986950000000 pi 

- 22384370986950000000 p2 
+ 49070938130697pip2 

- 67626498450000000000000000000 pi Ce (19) 

- 67626498450000000000000000000 pa Ce 
+ 163214120034000000000000pip2 Ce 

+ 45346707000000000000000000000000 pip2 Cf . 

Further analysis shows that Ce = results in fc = 0, and the 
solution of the linear (with respect to Ce) equation implies the 
negative value of Cq according to ( fTSl l. Thus the two solutions, 
whose existence we claimed above, have their origin in the 
quadratic equation ( fT9] i. 

V. What If equations become too complicated? An 
Outlook 

Symbolic methods always suffer from one major deficiency 
- the "combinatorial explosion" of terms - which most of the 
time results in both, a massive, nearly insatiable demand for 
space, and an extraordinary high amount of computing time. 
This observation also applies to all Grobner basis methods: It 
is known that by necessity the problem of finding a Grobner 
basis in the general case needs an exponential amount of 
spac^lH. On the other hand, in recent years there have been 
major improvements with regard to the run-time behaviour of 
the original Buchberger algorithnQ which have resulted in 
faster and thus more powerful algorithms to compute Grobner 
bases. 

As we could see in the above examples, even moderately 
sized small circuits already give rise to mid-sized equations 1 
which again - by way of the methods from computer algebra - 
result in even bigger symbolic expressions for the sought-after 

"For more details and/or a mathematically more precise formulation of 
this statement, see |7|. 

'^See e.g. the results by J.-C. Faugere in |8| and 19]. 
Indeed we must not regard the results, which have been presented up until 
now, as large - the maximal computing time used in the above examples in 
the worst case, which was the calculation of all the solutions of the Chebyshev 
RC filter, was ten minutes. 
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design parameters. Thus the "new toolbox" computer algebra, 
the usage of which I try to recommend with this article, seems 
to be doomed from the start. 

Still, bad behaviour in general does not mean that the in- 
struments from computer algebra always behave badly - there 
might be classes of problems for which results are computed 
fast, or can be written down very economically. As electrical 
engineers know, there are close (but in the author's opinion by 
far not completely understood) connections between a circuit's 
topology and the algebraic structure of its network function, 
which may be taken advantage of with regard to design tasks. 

Take on the other hand researchers in symbolic methods, 
who are accustomed with the phenomenon of swelling ex- 
pressions arising in the analysis of circuits. As a workaround 
a number of them have proposed so-called "sequences of 
expressions It is interesting to note that even electrical 
engineers without experience in methods of symbolic analysis 
may already be well versed in this "sequential approach", if 
only they knew classic filter theory, in particular the Cauer 
canonical forms of LC-circuits. 

In form of an example, let us show, what we mean. Suppose 
we are given an LC-ladder network, consisting of three 
inductances and three capacitances, as shown in figure 5. 



L1 



CI ; 



L3 



Fig. 5. LC ladder in Cauer canonical form 

As is well known its impedance Z{s) can be described by 
a continued fraction 



Z{s)^Lis + 



1 



(20) 



Cis 



L2 s 



Cos 



L?. s 



That, in fact, is just short hand for a sequence of expressions. 
This can be seen by deducing ( l20l i with a process of alterna- 
tively combining impedances in series and in parallel. Thus 
we have 

Z(s) =iis + Zi(s), 



where Zi (s) denotes the impedance of the ladder circuit, with 
Li having been removed. Further: 

1 _ 1 1 1 



(21) 



^See (To). 



where Z2{s) is the impedance of the circuit, that is left, when 
Li and Ci have been removed. Continuing like this, taking 
away one more circuit element in each step, this process can 
be repeated, until we get to one last impedance Z^{s) (which 
would just consist of C3). 

Reinserting the various Zi we get 



Z{s) 



Lis + Zi{s) 
Li s 



Cis- 



and thus are able to calculate Z{s). 

When conversely one wants to realize an LC trans- 
impedance function as an LC ladder, classic theor}0 tells us to 
expand the given rational function into a continued fraction (an 
algorithm which in case of rational functions always comes to 
an end after a finite number of steps) and read off the element 
values directly from this expressiorF^ which according to 
theorems in classic filter theory all are positive. 

If we want to bring together the capabilities of a CAS for 
solving algebraic equations and the sequence of expression 
approach in our example, we first observe that we can easily 
calculate the form of the intermediate impedances Zi, i.e. we 
have 



1 



7^ = —, Ziis) = s + Z^is) 
C3 s *s 



1 



C2S + 



1 



Z^is) 
- a2S^ 



1 



ass 



5 -I- ass^ -I- ais' 



and finally 



Z{s) = Lis + Zi{s) 

_ (as • Li)s^ + (03 • Li + 04)3^ + (ai • Li + 03)5^ + 1 



ais 



with * and the denoting polynomials in the circuit para- 
meters Ci , C2 , C3 , L2 , L3 . 

A sensible design task would now be to specify an 
impedance function 



Zdes{s) 



k {AqS^ + AjS^ + A2S^ + I) 

A5S^ + Ass^ + Ais 



(22) 



with real numbers k,Ai,...,AQ, and to demand Z{s) — 
Zdes{s). Analogously to the above example of the BIT am- 
plifier, by matching coefficients, this results in six nonlinear 
equations, which have to be solved for the variables ai, . . . , 05 
and Li. Here we only present the (easy) solution (determined 
by Mathematica again) for the first step, which is 

Ai-Ae ^ As- Ae 



a2 = A2 



05 = 



k ' 

As ■ fc 
A. 



A. 



fli 



04 = A4^ 
Ax _ A-i 



A. 



•'See e.g. fTTI . 

'*Again, nowadays a computer algebra system might be of help for 
calculating this expansion without much effort, see 1121 . 
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To calculate further element values, this procedure needs to 
be repeated for the other Zi. Thus one has to solve next the 
equations resulting from Zx{s) = Ci ^ + zl{s) ' where (s) 



should still be described in terms of the coefficients ai 
Zi{s) is given by the ansatz 



64 5^ + 62 5' + 1 

63 + fei s 



and 



(23) 



with new unknowns 61, . . . , 64. This has to be continued until 
we have equations for each set of coefficients of Zi in terms 
of the coefficients appearing in plus additional equa- 

tions for the element values. Starting from those equations, 
successive reinserting leads to the desired expressions for the 
element parameters. 

One has to concede that - compared to the traditional 
continued fraction expansion - this procedure looks very 
cumbersome. But in fact, if one looks a little bit closer, it 
is equivalent to repeatedly doing polynomial division as part 
of a continued fraction expansion. Furthermore our approach 
might be generalized to other impedance functions which 
need not come from ladders, but still ultimately result from a 
sequence of series, respectively parallel reductions applied to 
the initially given circuit, i.e. the process of replacing two 
impedances in series, or in parallel by a single equivalent 
impedanc^lZl Beyond this there is a further generalization to 
handle Delta-Wye and Wye-Delta transformations. 

On the other hand by results from Epifanov and Truemper, it 
is known that every planar two terminal graph can be reduced 
to a single edge by these four kinds of transformations (plus 
removing loops and pendant edges). Moreover there is an 
algorithm, developed by Feo and ProvarTl which calculates 
the sequence of necessary graph transformations. Since to 
each of these transformations there corresponds an equation, 
which gives the value of the resulting element(s) in terms 
of the values of previously existing elements, there always 
is a sequence of equations which leads from the initially given 
element values to the driving point impedance of a planar two 
terminal circui 

Circuit design for (midsized to large) planar two terminal 
circuits, which consist of impedances only, thus becomes the 
task of "reversing" these equations, in an analogous manner 
to what we did with LC ladders. 

The author of this paper firmly believes, that computer al- 
gebra systems will be a valuable tool in tackling this particular 
problem. The challenge remains to find classes of circuits for 
which the sequence of reversed equations can be effectively 
calculated, or to prove - what might be also possible, and 
would be an interesting result in itself - that, besides the 
classic, continued- or partial-fraction-topologies attached to 
the names of Cauer and Foster, no such further classes exist. 

'^However we need to be careful with the claim for generality, since if we 
combine e.g. two resistances i?i , R2 in series to a single resistance R, from 
the resulting equation R = Ri + R2 , when only given the value of R, we 
cannot deduce uniquely the values of the Ri. 

'^As a reference for these results in graph theory, see the original paper 

(H) 

"This algorithm has been successfully implemented in a student's minor 
thesis 1141 at the TU Braunschweig in 2003. 
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